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ABSTRACT 

We investigate the electromagnetic (EM) counterpart of gravitational waves (GWs) 
emitted by a supermassive black hole binary (SMBHB) through the viscous dissipation 
of the GW energy in an accretion disk and stars surrounding the SMBHB. We account 
for the suppression of the heating rate if the forcing period is shorter than the turnover 
time of the largest turbulent eddies. We find that the viscous heating luminosity in 
O.IMq stars can be significantly higher than their intrinsic luminosity, but still too low 
to be detected for extragalactic sources. The relative brightening is small for accretion 
disks. 
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1 INTRODUCTION 

The coalescence of supermassive black hole binaries (SMB- 
HBs) generates gravitational waves (GW) which are a pri- 
mary source for the proposed Laser Interferometric Space 
Antenna (LISAq). SMBHBs are inevitable outcomes of 
galaxy mergers. Spatially-resolved active galactic nucle i have 
been observed ([K omossa ct al. 2003; Bianchi ct al. 
iGreen et all bOloTlKoss et all I2OIII : iFabbiano et alT 



In addition, spectroscopic surveys (jComerford et al. 




ISmith et aU I2OI0I : iLiu et all l2010bl ') and observations that 
combine ground-based imaging show numerous systems con- 
taining compelling SMBHB candidates with pc t o kpc sepa- 
ratio n s ((Rodriguez et alJl2006l:[Liu et al.ll2010al : IShen et all 
I2OIII : iFu et al.ll201ll : iMcGurk et all bOlll l.llvdrodvnamic 
simulations of galaxy mergers also predict SMBHB pair for- 
mation (e.g.LE scala ct al.'2004, 2005; Di Matteo ct al. '2005"; 
Robertson et al. 2006; Hopkins et al. 2006; Callegari et al. 



20091 : IColpi et al.ll2009l : TBlecha et al.ll2012l ). 

Electromagnetic (EM) counterparts to GW sources 
complements the GW detection by determining the host 

J alaxy redsh i ft and the environ ment of the sources 
Kocsis et all I2OO6I : IPhinnevI bOOgl ). A large variety of 
EM signatures have bee n proposed to accompany the co- 
alesc ence of SMBHBs ISchnittmanl I2OIII : iHaiman et al.l 
I2OO9I ). In the pre-merger phase, the torques of the 
SMBHB excavates a hollow region in the disk and 
leads to p eriodic accret i on across the gap on the orbital 
timescale JCuadra et al. _ 2009 : iMacFadven fc Milosavlievig 
l2008l : iHavasaki et alj|2008l ') .After the merger, the recoil of 
the black hole remnant and its sudden mass loss due to the 
final GW burst produ ce shocks in the accreti o n disk which 
lead to EM signals (|Bode fc PhinnevI l2007l : iLippai et al] 



http://lisa.nasa.gov/ 



20081: ISchnittman fc Kroli^ I2OO8I: IShields fc BomiiiiS I2OO8I: 



O'Neill et al.ll2009l : lRossi et al.ll2010t ). The recoil of the black 
hole remnant changes the tidal disruption rate of stars due to 
the refilli ng of the loss cone an d the wandering of black hole 
remnant ( JStone fc Loebll2011al lbl: lLi et aLlbOl g). Finally, the 
infall of g as onto the black hole remnan t produces an EM 
after glow (|Milosavlievic fc Phinnevll2005l : iTanaka fc Menoul 
|2010D . 

In this paper, we consider the viscous dissipation of 
GWs generated by a SMBHB in a neighboring gaseous 
medium. In particular, the velocity shear induced by GWs 
in the gas is damped by viscosity. The dissipated GW en- 
ergy turns into heat, and produces an electromagnetic flare. 
Unlike other EM counterparts, the brightening here follows 
promptly withi n a few hours to day s after the coalescence 
of the SMBHB (JKocsis fc Loebll2008l '). The effect provides a 
unique test of general relativity for the interaction of GWs 
with matter. In § [5] and [3] we investigate GW dissipation 
in a gaseous accretion disk and stars in the vicinity of the 
SMBHB. We examine the suppression of the effect if the 
forcing period is shor ter than the turnover time of the largest 
eddies (|Krolikl I2OI0I V in an alogy to a similar treatment of 
tidal heating in binary stars (l Zahnll966l : lGoldreich fc KeelevI 
[1977). Finally, we discuss our conclusions and their implica- 
tions in § U 



2 METHOD 

We start by presenting our approach for estimating the GW 
heating inside an accret ion disk and stars due to turbulent 
viscosity. Following Kocsis fc Loebl (|2008l ). we approximate 
the GW luminosity by matching the Newtonian inspiral lu- 
minosity prior to merger (i < 0), the peak luminosity at the 
merger (t — 0) and the decay luminosity afterwards (i > ti). 
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where ti can be fixed from this matching procedure. Specif- 
ically, in the Newtonian inspiral regime, the luminosity is 

^ 32 G'-'mV^ ,,^ 

ivGWinspiral - y-^— ^5— , (1) 

where M = Mi +M2 is the sum of the masses of the SMBHB 
members, /x = M1M2/M is the reduced mass of the SMBHB 
and a is the separation between the SMBHB, which can be 
expressed as 

-,1/4 



256 G^ 2,, 



t) 



(2) 



assuming a circular orbit. The peak l uminosity i s ap- 
proximated from num erical simulations (JBerti et al.l 120071 : 
iBuonanno et al.|[2007l ') as 



and the ringdown luminosity is set to be 



(3) 



r r ( c{t~ti)\ .^. 

^GW ringdown — ^GWpoakexpl — — I , (4) 

where Rg = GM/c? is the gravitational radius of the 
SMB HB. The peak luminosity is modified by a factor of 
two (JBerti et al.ll2007l : [Buonanno et al.l 120071 ') due to differ- 
ent magnitudes and orientation of the spin of the SMBHB. 
In this paper, we assume the masses of the two black holes 
are the same. 

With the approximated expression of GW luminosity 
as a function of time, the dissipation of GW energy inside a 
viscous medium can be calculated by solving the weak-field 
Einstein equation (Hawking 1966; Weinbe rg,, 19721 ): 



16-nGr] 

Cheat = 7, ecw, 



(5) 



where fiheat is the dissipation rate, r) is the dynamical vis- 
cosity and ecw is the GW energy density, ecw can be 
obtained from ecw ~ y(^) ^°^ , where 9 is the angle 
relative to the total angular momentum vector, Y[6) = 
5/2[sin*(6l/2)-|-cos*(6l/2)]. We use the average value (Y) = 1 
below. With Lew derived, the only unknown parameter is 
the dynamical viscosity of the medium that the GW passes 
through. The dissipation rate of the GW energy gives the 
heating rate of any gaseous medium such as an accretion 
disk and stars. 

Next, we estimate the dynamical viscosity for stars. We 
use stellar models produced by Modules for Exper iments 
in Stellar Astrophysics (MESAQ) (|Paxton et al.ll201ll '). a ID 
stellar evolution code, and we consider stellar models, whose 
properties are included in Table [l] We associate the dynam- 
ical viscosity with the mixing length theory diffusion coeffi- 
cient, which is directly provided in the simulated models by 
MESA. When the period of the driving force is smaller than 
the largest eddy turnover time, the eddy viscosity depends 
on the ratio of the period to the largest eddy turnover time 
in one of two possible ways: 



^.^,,nin[(^),l] 



(6) 



■q = rii mm 



(^-^\\l\ 



(7) 



LV27rr,/ 

where rji is the intrinsic viscosity in the absence of shear force 
with short period, rj is the largest eddy turnover timescale 
and TGw is the shear force period, which is calculated as 
2-k/ujgw, where ojgw = 2^J GM j a? in the inspiral phase a < 
6-Rg and 0. 25/(6* M/c^)] after the ringdown, and extrapolate 
linear ly during the transition according to iBuonanno et al.l 
1120071 ) . The viscosity scaling given by Eq. ([§} is discussed in 
ZahnI (ligedfigsgl: Zahn fc Bouchet (1989) and Eg . 



Goldreich fc KeelevI (|l977l ): iGoldreich fc NicholsonI l|l9& 



Observations are more consistent with Zahn's scaling for 
pulsating star s in the red edge of the instability strip 
(Gonczi 1982) , for tidal circularization of binary stars 
(Verbunt fc Phinnevlll995l : lMeibom fc Mathieull2005i ). while 
the damping of the solar p-mode oscillations is more consis- 
tent with the Goldrei ch's scaling (IGoldreich fc Kumarlll988l : 
IGoldreich erai]|l994 ). Recentlv. lPenev et al.l (|2009l ) studied 
turbulent visco sity in low mass stars using the perturba- 
tive approach of lGoodman fc OhI (| 19971 ) , taking into account 
compressible fluid and anisotropic viscos ity. Their simu- 
lation suggests a linear scaling. However, lOgilvie fc Lesurl 
(2012) found results more consistent with Goldreich's scal- 
ing when studying the limit of a low amplitude short oscil- 
lation period shear. We considered both scalings for stars in 
this paper. 

With the viscosity for stars and Lawit) in hand, the 
GW heating rate can be estimated using Eq. ((5)|. The EM 
luminosity increase can be estimated by solving the radiative 
transfer equation: 



t,(r)-A/(r) + A/(r) 



Cheat , 



Lg\ 



A/(r) dV, 



(8) 
(9) 



where A/(r) is the excess EM signal produced per unit vol- 
ume as a function of location in the star, Lgwh is the excess 
EM luminosity associated to GW heating, and tc{r) is the 
cooling time as a function of the location, which character- 
izes the time it takes for heat to travel to the surface. We 
estimate the latter by taking the integral of the minimum 
of the photon diffusion time, dr/c x [r{r) — (i?* — r) ^J" ], 
and the turbulent convection time, dr/vc{r), in each spher- 
ical shell inside the star, where the optical depth, T{r), and 
the convective velocity, Vc{r) are obtained from the MESA 
simulation, and Rt, is the radius of the star. 

Finally, we estimate the heating in accretion disks. We 
adopt the geometrically thin, optically thick, standard ac- 
cretion disk model, where the angular momentum transport 
is a ssociated with the int e rnal stresses due to turb ulence 
(Sha kura fc SunvaevI [l973l : iNovikov fc Thornelll973l ). Heat 
is dissipated locally by turbulent viscosity, and transported 
vertically outward by photon diffusion or advection. Specif- 
ically, the viscosity of the accretion disk is 

2 aP{r) 

3 n{r) ' 



Viir) 



(10) 



http://mesa.sourceforge.net/ 



where 0,'^{r) = GM/r'' is the a ngular velocity, a is a con- 
stant which we assume to be 0.3 (JKing et al.ll2007l ). and P is 
the total (gas-|-radiation) pressure in the a disk model, and 
gas pressure in the /? model. In these models, the physical 
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characteristics of the disk is fixed by the following parame- 
ters: the accretion rate in Eddington units ( rh), the radiation 
efficiency (e), and the SMBHB mass (M) (|Goodmanll2003l : 
iGoodman fc T an"2004Y We set m to be 0.1, e to be 0.1, and 
discuss the effects caused by different SMBHB masses. 

Similarly to stars, we account for the frequency depen- 
dence of viscosity when the period of the driving force is 
smaller than the largest eddy turnover time, and estimate 
the effective vis cosity according to the p erturbative methods 
as discussed in lGoodman fc OhI ((19971). Specifically, for in- 
compressible fluid with isotropic viscosity, the viscosity as a 
function of driving force frequency can be expressed in terms 
of the frequency spectrum of the average kinetic energy per 
unit mass. For accretion disks, where the Kolmogorov scal- 
ings may not be applicable to obtain the energy spectrum, 
we adopt the energy spectrum fr om recent magnetoh ydrody- 
namic (MHD) disk simulations. iFlock et al.l (|201ll ) present 
a full 27r three dimensional simulation on a stratified accre- 
tion disk, where the turbulence is driven by magnetorota- 
tional instability (MR I) and the kinetic spectra is obtained 
in the 4> direction, and lFromangj (|2010|) investiga t e the MRI 
in a shearing box with zero net flux. iFlock et al.l (|201ll ) and 
iFroma ng (2010) estimate the kinetic energy spectrum expo- 
nent to be 11/9 and 1.5, respectively. We estimate the vis - 
cous heating in accretion disk with the iFlock et al.l (|201lh 
exponent as well as the Kolmogorov scaling exponent 2. 

Similarly to the calculation for stars, the GW heating 
rate can be estimated using Eq. ((5]). The corresponding EM 
signals can be est imated by solving the radiative transfer 
equation following iKocsis fc Loebl l|2008l ): 



tc(r) ^ AF(r, t) + AF{r, t) = Hehoat (r, t), 
at 



iGWH(i) = 



2TvrAF{r,t) dr, 



(11) 
(12) 



where F{r) is the excess EM flux due to GW heating in 
the accretion disk, H is the scaleheight, Lgwh is the cor- 
responding excess EM luminosity, and tc{r) is the cooling 
time. Here, we assume the disk is face on, and account for 
the different light-travel time from different annuli in the 
disk. The brightening can be somewhat larger in an inclined 
or edge-on configuration (by up to a factor of ~ 3) where 
the peak GW flux is observed coinc identally at the inner and 
outer radii along the line of sight l|Kocsis fc Loebll2008l l. 



3 RESULTS 

First, we consider the GW heating of nearby stars. As an ex- 
ample, we examine the GW heating light curve for a 0.1 Mq 
star (stellar model 2) surrounding an M = 10^ or an W^Mq 
SMBHB, respectively. Using Eqs. ((8]) and ([9]), we calculate 
F{t) and plot the GW heating light curve in Figure [T] We 
assume that the star is located at d = 5 tidal radii from 
the SMBHB (corresponds to 320 and 15 Rg for a lO'^ and a 
10^ Mq SMBHB, respectively). Note that since the GW lu- 
minosity is proportional to {d/Rg)~'^, the GW heating effect 
is much larger around more massive SMBHBs because the 
viscosity suppression for a high mass SMBHB is smaller. 

Figure [T] shows that the excess luminosity of the star 
surrounding the 10^ Mq SMBHB is much higher than the in- 
trinsic luminosity of the star {L 
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Figure 1. The light curve of a GW heated star (based on star 
model 2 with an intrinsic luminosity: L = 2.6 X 10'^'' erg s~^). The 
time axis is in units of Rg/c, and is shown on a logarithmic scale 
at both negative and positive values (causing the discontinuity at 
t = 0). The star is located 5 tidal radii away from the SMBHB 
(320 and 15 Rg for a lO'' Mq and 10^ Mq SMBHB, respectively). 
The black line indicates the GW luminosity scaled down by 25 
order of magnitude in order to fit in this figure. The red and 
blue lines indicate the light curve of a star surrounding a 10 
and a 10® Mq SMBHB, respectively, with solid and dashed lines 
corresponding to the viscosity dependence with (tgw/2t;) and 
(TQw/27rT;)^, respectively. The light curve closely tracks the GW 
light curve. Interestingly, the peak luminosity surrounding the 
10® Mq SMBHB is much higher than the intrinsic luminosity of 
this star. 



the net dissipated GW energy can exceed the gravitational 
binding energy near the stellar surface, and could generate 
a stellar wind. However, as the viscosity is strongly sup- 
pressed in the stellar interior (^^^ "C 1 for r < 0.99Ratar), 
the heating effect is negligible to the star as a whole. In 
addition, these stars are very faint; the absolute peak GW 
heating luminosity in the star is typically too faint to be 
observed outside of the Galaxy. 

Since the turnover time of turbulent eddies is much 
longer in the interior of the star than that at the surface, the 
energy is mostly dissipated at the surface. Since the cooling 
time near the surface (~ 200 s) is short compared to the peak 
GW timescale (~ WRg/c ~ 500 s A/bh/IO^ Afg), the light 
curve of the star closely tracks the luminosity curve of the 
GW. When the GW driving period is shorter than the eddy 
turnover time, the viscosity caused by the eddy depends on 
the ratio tqw/ti, where the exact scaling is uncertain as 
discussed in § [2l For stars surrounding a 10® Af© SMBHB, 
the differences between the two scalings are smaller as the 
period of peak GW emission for this SMBHB mass is more 
comparable to the surface eddy turnover time in a 0.1 Mq 
star. 

To examine the influence of the GW heating in different 
types of stars, we consider stellar models of different stellar 
masses and ages as included in Table 1. We include the ex- 
treme CEises with 0.1 Mq and 100 Mq stars. We plot the ratio 
of the peak heating luminosity to the intrinsic luminosity for 
different stellar models in Figure 2 with Zahn's scaling. We 
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Table 1. Properties of stellar models 



No. Mass Metallicity Radius Luminosity Age 

(Mo) (Z) (i?,0) (Lq) (yrs) 



1 


0.1 


0.16 


2 


0.1 


0.16 


3 


0.1 


0.16 


4 


0.1 


0.01 


5 


0.1 


0.01 


6 


0.1 


0.01 


7 


100 


0.04 


8 


100 


0.04 


9 


100 


0.04 



3.3 


0.00079 


2 X lO'' 


3 


0.00066 


5 X 10*5 


0.57 


5.5 X 10-^ 


2 X 10^ 


2.4 


0.43 


2 X 10-* 


0.44 


0.022 


5 X 10** 


0.12 


0.0012 


2 X 10^ 


21 


1.4 X 10^ 


1 X 10-* 


36 


1.7 X 10*^ 


1 X 10*5 


960 


2.1 X 10^ 


2 X 10** 




10" 10- 10'" 10" 10'" 10" 10'" 

Mbh(M0) MBH(Me) Mbh(Mq) 



Figure 2. Ratio of the peak GW heating luminosity to the in- 
trinsic stellar luminosity. The horizontal axis shows the mass of 
the SMBHB, and the vertical axis plots the distance (d) between 
the star and the SMBHB in units of the tidal radius (rt). First 
row: model 1, 2, 3; second row; model 4, 5, 6, third row: model 
7, 8, 9. Solid black line indicates where the distance between the 
star and black hole binary is 6 Rg, the radius of the innermost 
stable circular orbit (ISCO) around a non-spinning black hole. In 
the last panel, the points in the figure lie out of 6 i?,g, and so the 
black line is not shown. The first two rows correspond to 0.1 Mq 
stars with metallicity Z = 0.16 and Z = 0.01 respectively, and 
the last row corresponds to 100 Mq stars. GW heating is most 
significant for high metallicity low mass stars. 



find that the influence of the GW iieating is more significant 
as the metallicity of the star increases, and GW heating is 
not significant for very massive (M* > IOOMq) stars. 

Next, we discuss the heating effects in accretion disks. 
For a and /3 disks, we solve Eqs. (|lip and (|12|) for the heat- 
ing flux, and plot the heating light curve of the disk due to 
GW heating in Figure 3. The accretion disk is punctured 
with an inner hole. This geometry is essentially "frozen" 
during the final GW merger timescale with a gap r adius 
> lOOM for a-disks dMilosavlievic fc Phinnevll2005l ). Re- 
cent MHD simulations bv lNoble et al.l (|2012l ) indicate that 
the stresses may be enhanced in a binary, such that gap de- 
coupling occurs further in, at 20Rg. We optimistically adopt 
this value for our estimates, which implies a larger heating 
rate than that for a larger gap radius. We integrate over the 



10 



10 



.a 
^10 



10 



10 



5 


y 


/^II— -^ 




/^^^^^^^— -— ^^^ 


6 




-'•■'"■ \\ 


7 




Flock tc<< r/c 

Flock /i disk 

Kolmogorov tc << r/c 

Kolmogorov a disk 

Kolmogorov disk 


1 v\ \\ 


8 




, 


\ \\ \\ 



-10000 



-100 




time (Rg/c) 



100 



10000 



Figure 3. The excess luminosity relative to the disk luminosity 
due to GW heating on an accretion disk (inner disk truncated at 
20 -Rg) before (t < 0) and after (t > 0) the binary coalescence 
event. The time axis is shown on a logarithmic scale at both 
negative and positive values (in units of Rg). The SMBHB mass 
is W'^Mq. Solid lines corresponds to the frequency dependence 
{tq\j\/ /2iTTi)~9' derived according to the energy spectrum of accre- 
tion disk based on MHD simulations by Flock et al. (2011), and 
the dashed lines correspond to the scaling (rQw/27rr;)^, assuming 
Kolmogorov turbulence. 

accretion disks between the inner and outer boundary. We 
set the latter to 2 x W^Rg, but this value does not infiu- 
ence our result as the heating in the outer accretion disk 
is negligible. We include the different light travel time from 
different accretion disk surface elements along the line of 
sight. Our calculation of the heating in the accretio n disks 
improves the simplified treatment of iKocsis fc Loebl (|2008h 
by including the dependence of viscosity on the ratio of the 
GW driving period to the largest eddy turnover time, which 
suppresses the dissipation of GWs. We consider two cases in 
th is plot. Following the p erturbative turbulence derivation 
bv lGoodman fc OhI ([1997|), the power-law index is 2 for Kol- 
mogorov turb ulent scaling , and ^ according to MHD disk 
simulation bv iFlock et al.l (|201ll ). The eddy turnover time 
increases rapidly as the radius increases, and so the suppres- 
sion of the GW heating is less significant for disks truncated 
closer to the SMBHB. Therefore, the heating luminosity is 
more significant for disks that are truncated closer to the 
SMBHB. 



4 DISCUSSION 

In this paper, we considered the dissipation of GWs in an 
accretion disk or stars surrounding a SMBHB. We have 
found that the GW heating luminosity of the accretion disk 
and stars are low, and make no significant EM fiare rela- 
tive to their intrinsic luminosity except for low mass stars 
(~ O.IMq). The integrated excess luminosity from heated 
low mass stars is too low to be observed in galactic nuclei 
as they are faint. Assuming a Bahcall-Wolf distribution of 
stars or assuming a collision timescale larger than 1 Myr, we 
find that only a few stars are expected to be within 5 tidal 
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radii of a coalescing SMBHB, where the GW heating effect 
is significant. Therefore the overall brightening of the stellar 
cluster is negligible. 

In order to be heated significantly by GWs, the stars 
need to be close to the SMBIfB. One possible avenue is 
that stars get caught in mean motion resonances (such as 
Trojan re sonances) and mov e inwards as the SMBHB merge 
l|Seto fc M uto 20f0l: ISchnittm an 20i0). This is only effective 
for SMBHB with an unequal mass ratio q < fO~^; the stars 
get ejected before the coalescence otherwise. Another possi- 
bility is for stars to get captured or form in the outer parts 
of accretion disks, and migrat e inwards by pro cesses a nal- 
ogou s to planetary migration ([Miralda-Escude fc KoU meierl 
I2OO5I : iKaras fc Subrll200ll : lLevinll2007h . 

We assumed that GW energy is dissipated locally 
through turbulent viscosity. The damping of shear stress by 
eddy viscosity in stars was found to be consistent with ob- 
serva tions in the context of the tidal circularization of bina - 
ries l|Verbunt fc PhinnevI [l995l : iMeibom fc Mathieull2005h . 
The underlying accretion disk model is uncertain since the 
disk structure is u nstable to both the rmal and viscous insta- 
bilities. Recently, Blacs ct alj (|201l|) found that radiation- 
dominated disks differ significantly from the standard disk 
models, where the dissipation associated with the turbulent 
CEiscade and radiative damping dissipate energy non-locally. 
It remains to be seen whether the GW heating effect is more 
prominent in alternative disk models. 
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